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Csj ; Abstract 

In analogy with ordinary BCS superconductivity, we identify possible pairing instabilities 

\^ • with poles in a certain class of Green functions of the two-dimensional repulsive Hubbard 

model, which normally signify bound states. Relating the gap to the location of these 

D I poles, we find that these instabilities exist in the range of hole doping between and 0.24. 

'^ ' The magnitude of the gap can be calculated without introducing an explicit cut-off, and 

j^ . for reasonably large coupling U, the maximum gap is on the order of 0.01 in units of the 
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I. INTRODUCTION 

The microscopic physics underlying high Tc superconductivity in the cuprates is 
beheved to be purely electronic in origin, in contrast to ordinary superconductors 
where the attractive mechanism is due to phonons. Strongly correlated electron 
models such as the two-dimensional Hubbard model have been proposed to describe 
it. The Hubbard model simply describes electrons hopping on a square lattice subject 
to strong, local, coulombic repulsion. Since it is known that the condensed charge 
carriers have charge 2e, and thus some kind of Cooper pairing is involved, the main 
open problem is to identify the precise pairing mechanism. Several mechanisms were 



proposed early on, in particular mechanisms based on spin fiuctuations[lH3|, charge 
fluctuations|5|-l7|], and other more exotic ideas such as resonating valence bonds 8|. 
A more recent work studies the possibility of superconductivity at small coupling, 
where here the mechanism goes back to ideas of Kohn-Luttinger|9|, [lO|. 

Since the Mott-insulating anti-ferromagnetic phase at half-filling is well under- 
stood, much of the theoretical literature attempts to understand understand how 
doping "melts" the anti-ferromagnetic (AF) order, and how the resulting state can 
become superconducting. This has proven to be difficult to study, perhaps in part due 
to the fact that AF order is spatial, whereas superconducting order is in momentum 



space. (For a review and other refereces, see [llj.) Furthermore, superconductivity 
exists at reasonably low densities far from the AF order at half-filling, which suggests 
that that one can perhaps treat the model as a gas, with superconductivity arising 
as a condensation of Cooper pairs as in the BCS theory, and we will adopt this point 
of view in the present work. 

Since there is as yet no consensus on the precise pairing mechanism in the cuprates, 
it is worthwhile continuing the search for new ones. In this paper we take a very 
conservative approach within the Hubbard model, wherein we do not postulate any 



particular quasi-particle excitations that would play the role of the phonons, and we 
ignore the AF order at very low doping. In a field-theoretic approach to the BCS 
theory, the pairing instability arises from the sum of the ladder diagrams for phonon 
exchange, as shown in Figure [H The sum of these diagrams has a pole at energy 
equal to the gap, and this pole is sufficient to cause the instability, and signifies 



gapped charge 2 excitations. This will be reviewed in section III. (See e.g. 
> . > . > . > . > . > 



12 



13|). 
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FIG. 1: Feynman diagrams that lead to the BCS pairing instability, where dashed lines 
are phonons. 

Suppose that the electron-phonon interaction were treated as an effective, short- 
ranged electron-electron interaction, so that exchange of a single phonon were re- 
placed by an effective 4-electron vertex. Then the diagrams in Figure [H become the 
diagrams in Figure |2j In this work we explore the possibility of pairing instabili- 
ties arising from certain similar classes of Feynman diagrams in the Hubbard model. 
Whereas the diagrams in Figure [2] lead to the BCS pairing instability for bare at- 
tractive interactions (see below), for repulsive interactions they merely screen the 
strength of the Coulomb interaction. However, as we will show, the sum of another 
class of diagrams shown in Figure [3] do indeed exhibit poles possibly signifying pair- 
ing instabilities very near the Fermi surface. An obvious criticism of the present work 
is that for the cuprates, the coupling is large, and one should not trust perturbation 
theory. In answer to this, the diagrams we will focus upon can be summed up, and 
have a well-defined limit as the coupling goes to infinity. Furthermore, since we are 
focussing on sums of diagrams that lead to poles in the effective pairing potential, it 



is possible that these singular diagrams dominate the perturbative expansion. Cer- 
tainly our calculation can be improved upon; we wish here to present a possible 
pairing mechanism in its simplest form. For instance, we ignore the effect of interac- 
tions on the quasi-particle energies, and take ^k to be that of the free theory, ignoring 
self-energy corrections. 

The remainder of this paper is organized as follows. In the next section we define 
the Cooper 'pairing' Green functions of interest. In section III we first review how 
the Cooper pairing instability manifests itself as poles in such Green's functions in 
the BCS theory. The remainder of this section specializes to the Hubbard model, 
where analogous poles are found for a different class of diagrams. The gap is related 
to the location of these poles, and solutions are found numerically. We find non-zero 
gap solutions in the anti-nodal directions in the range of hole doping < h < 0.24 
with a dome-like shape. The maximum of the gap is approximately A/t = 0.01 and 
occurs around h = 0.18. 
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FIG. 2: Feynman diagrams contributing to V^jj 
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FIG. 3: Feynman diagrams contributing to V^g' 



II. GREEN FUNCTIONS AND THE EFFECTIVE COOPER PAIR PO- 
TENTIAL 

We study the two-dimensional Hubbard model, with an additional next-to-nearest 
neighbor hopping t', with hamiltonian: 

H = -t ^ cl,aCj>-^' Y^ c{aCj,a + uYni^^ni^^ (1) 

<i,j>,a="\,i <i,j>',a i 

where i,j label sites of a two-dimensional square lattice, < i,j > and < i,j >' denote 
nearest and next-nearest neighbors respectively, and rii^a = cj^Ci^a- Our convention 



t^ _ J^ I J, 



IS 



for neighboring interactions is e.g. ^^i jyclcj = c[c2 + CgCi + such that H i 

hermit ian. 

The lattice operators Cr-,^ where Fj is a lattice site, will correspond to the contin- 
uum fields 

M^,t)= /'^Ck,„(t)e^''-^ (2) 

J 27r 

The free hopping terms are diagonal in momentum space: 

Hiree = / rf^kWk ^ 4,aCk,a (3) 

with the 1-particle energies 

Wk = — 2t(cos(fc^a) + cos(A;j^a)) — it' cos(A;^.a) cos(A;j^a), (4) 

where a is the lattice spacing. We henceforth scale out the dependence on t and a 
such that all energies, temperatures and chemical potentials are in units of t. The 
interactions then depend on the dimensionless coupling g = U/t, which is positive 
for repulsive interactions. Henceforth C,]i = uj]f_ — fi is the 1-particle energy measured 
relative to the Fermi surface, where // is the chemical potential. For the most part, 
our analysis will be for t' = 0, unless otherwise stated. 



In order to probe possible Cooper pairing instabilities, we consider the Green 
function {^I{xi)^l{x2)4'i{x3)^^{x4)) where x = (x, t). The Fourier transform in 
both space and time of these functions are correlation functions of the operators 
Ck,«(€) — / ~^^~*^*'^k,a(^)' ^^^ their hermitian conjugates. We study Green func- 
tions specialized to Cooper pairs: (cl,^,|(f)c[.,| (c')c_k4,(C)Ckt(?))- It should be em- 
phasized that since the above Green functions are just Fourier transforms of the 
spatial/temporal Green functions, although the C, are energy variables, they are not 
necessarily 'on-shell', i.e. ^ is not necessarily ^^ (borrowing the relativistic terminol- 
ogy). The only constraint is energy conservation: ^ + ^ = ^' + ^'. 

There are two important types of quantum corrections to the vertex, which we 
denote as V^g and V^g , where, for reasons explained below, sc refers to screened 
and ex to exchange. V^^ is defined as 

Vir^ = (c^-k't(?)4,(r)c-H(ackt(«))llc (5) 

where trunc refers to the truncated Green function, i.e. stripped of external propaga- 
tors and energy-momentum conserving delta functions, and sc refers to the diagrams 
in Figure [2l Incoming (outgoing) arrows correspond to annihilation (creation) oper- 
ator fields. The other class of diagrams are shown in Figure [3] and define V^g^ as in 
eqn. ([5]). 

When 'on-sheir, i.e. .^ = .^ = .^k and C,' = C,' = ^k', then these truncated Green 
functions contribute to the matrix element, i.e. form-factor, of the integrated inter- 
action hamiltonian density: 



V{k, k') = j d2x(-k' t, k' ; |Hint(x) |k t, -k i) 



(6) 



Since there is no integration over time in the above equation, ^k and ^k' are not 
necessarily equal. To lowest order, V = g. The form factor V^'^'^^'ik, k') corresponding 
to the diagrams of V^g is equal to V^^ (as computed below) with ^ — )■ (^k + 



.')/2Q, 



.^k')/2(l4l' whereas \/(°'') is simply V^g^ placed on shell. For both, one has the 
necessary symmetry: ^(k, k') = V(k',k). 

The evaluation of V^^"^'*^^ at finite density and temperature is standard, however for 
completeness we provide some details. Consider first V^^^^ . These diagrams factorize 
into 1-loop integrals and form a geometric series. There is no fermionic minus sign 
coming from each loop since the arrows do not form a closed loop. Momentum 
conservation at each vertex gives a loop integral that is independent of k, k': 

where L^**^^ is the one-loop integral: 

where T is temperature, f„ is a fermionic Matsubara frequency, Un = 27r(n + |)T 
with n an integer, and ^p = Wp — /i where fi is the chemical potential. One needs the 
following identity: 



where /(.^) = l/{e^''^ + 1) is the fermionic occupation number, and we have used 
,^p = ,^_p. The above identity is valid before analytic continuation from imaginary 
to real time, i.e. when 2^ is twice a fermionic Matsubara frequency. The final result 
is then: 

^'"'K'4/|^(«3^)'-''(W2r) (10) 

where t] is infinitesimally small and positive. Integration is over the first Brillouin 
zone, —71 < p,j.y < IT. 

The diagrams for V^g^ , though simply an 'exchanged' version of those for V^^^^ , 
have a rather different and more complicated structure. Here there is a fermionic 



minus associated with each loop since the arrows form a closed loop. Momentum 
and energy conservation at each vertex now leads to a loop integral that depends on 
q = k' - k and Ag = f - ^: 



where 



Vr(q.A,) = ^3^ (11) 



^-(.-.)-E/|^(^)(: 



(27r)2 XiUn-^Y,) V^C + ^^n - Cp+q 

Now one needs the identity 



(12) 



"^ V^'^n - 4p/ V^^n + «t^m " ^p+q/ ^^m + ^p " 4p+q 

which is valid for Um twice a bosonic Matsubara frequency. After analytic continua- 
tion ium. — 7- A^ + i?7, one has 

L<"'(q. A,) = /|H_ ( /Kp)-/(W/ ) (14) 

i (27r)^ V^«+^p-^P+q + ^^/ 



In the above formulas it is implicit that V'-'^'^'"^-' is defined by the real part of L^^'^''^^^ 

,{ex 
eflf 



in the limit 77 — )■ 0. In the limit of zero density, i.e. / = 0, note that Vg^ = 



whereas V^g is non-zero. 



III. COOPER PAIRING INSTABILITIES AS POLES IN THE EFFECTIVE 
POTENTIAL 



It is well-known that the Cooper pairing instability of the BCS theory can be 
understood as a pole in V^^ 12|, [l3|. This can be seen by inserting a complete 



set of states between the pair creation and annihilation operators in eqn. (E]), and 
noting that poles in the energy ^ signify bound states of electric charge 2, which 



are interpreted as the Cooper pairs. To see that this gives the correct result for the 
gap, let us leave behind the Hubbard model and specialize the coupling g to that of 
the BCS theory, where it arises from the interaction with phonons. For simplicity, 
g = —\g\ is taken to be negative, signifying attractive interactions, in a narrow region 
near the Fermi surface |^| < lod-, where lob is the Debye frequency, otherwise zero. We 
approximate the integral over p near the Fermi surface as f p{^p)d^p ~ p(0) f d^p, 
where p{^) is the density of states. Letting C, = A/2, in the limit of zero temperature, 

L''-' - ^ f / dC.1^7^ - [ dCp—^] ^ p(0) log 



o<^p<i^D A/2 - ^p J-UD<^p<o A/2-,^py \2ljd 

(15) 

(The first integral should be interpreted as the principal value.) There is a pole 

in V^^ because both L^'^'^' and g are negative. The location of the pole is given 

by 1 + |5'|p(0) log(A/2ci;£)) = 0, which leads to the standard result for the gap: 

A = 2uj£) exp(— l/(|5f|p(0))). Note that there is a gap A for arbitrarily small negative 

coupling g. 

Returning to the positive g Hubbard model, i.e. with repulsive interactions, there 

are no poles in V^g since L^'^'^^ turns out to be negative at finite temperature and 

density. These quantum corrections merely screen the Coulomb potential. It is 

meaningful to define a screened coupling near the Fermi surface as g^ = V^^ {^ = 0). 

The extent of the screening is quite large, on the order of a factor of 10. For instance 

for g = 13, appropriate to the cuprates, g^ = 0.95 at half-filling. For the numerical 

analysis in this section, we fixed gji = 1. We wish to point out that at zero density, 

i.e. / = 0, and T = 0, L^'^^^ can become positive, which means that V ^"^ can 



be attractive. This was explored in 



I — I ' ' '^^ 

14l ]. and with the help of a hypothesized gap 

equation, one finds anisotropic solutions with properties suggestive of the pseudo- 
gap. 

We now turn to the investigation of the possibility of analogous pairing instabilities 

10 



in Vg^^ . Namely we look for poles in V^^'^ for positive g, which only exist if L^^^"* 
can be positive. For A^ = 0, the integrand in L^^^^ is always negative, thus L^'^^'^ is 
negative and there are no poles for positive g. However we are precisely interested in 
the case A^ ^ since A^ represents a difference of incoming and outgoing energies, 
and it is thus poles at non-zero values of A^ that could signify Cooper paired bound 
states. As we will see, there are indeed regions of the parameter space near the Fermi 
surface where L^'^^^ is positive. 

The integral defining L^'^^^ is rather delicate in comparison with L^'^^\ For q and 
Ag unrelated, the integral is easily evaluated numerically, and has a smooth, well- 
defined limit as ?7 — )■ 0. However physically we are interested in this function on-shell, 
i.e. the corresponding form factor ^(°^) with A^ = ^k' ~ ^k when q = k' — k. In 
this case there is a pole in the integrand when p = k. Thus the p integral should 
be understood as the Cauchy principal value (PV). Namely, inside an integral over 
a variable x, one has the identity: 

lim -^— = PV ( - ) - i7iS{x) (16) 

v^o+ X + irj \^ J 

Recall the PV is defined by excising a small region around the pole: J dp^ — )■ 
lim^^o (j^l~ dp^ + Jk^^s ^P^J ^^'^ similarly for / dpy. 

In order to better understand other features of the function L^*^^^ on-shell, it is 
instructive to study the one-dimensional version of L^^^\ appropriate to a chain of 
lattice sites, since here the integral can be performed analytically. This is described 
in the Appendix. (It should be emphasized that this exercise is not meant to capture 



the physics of the one-dimensional Hubbard model, which is exactly solvable (l 51] . but 
rather simply to gain intuition on the function L^'^^^ in two dimensions.) As shown 
analytically in the appendix, the on-shell L^^^^ actually diverges as the temperature 
goes to zero in one dimension. The reason is that in the limit of zero T, the occupation 
number / is a step function, and the f^ ^ dp piece of the PV is absent, leading it to 

11 



be ill-defined. The temperature T should thus be viewed as an infra-red regulator in 
one dimension. In two dimensions we will also not take T — )■ for analogous reasons. 
For arbitrarily small T, we verified that L^^^^ is well defined numerically as 5 — ?■ 0. 
As shown in the appendix, one can demonstrate analytically that there is a narrow 
region around the Fermi surface where L^*^^^ is positive and thus V^^g has poles. As 
we now show, the same is true in two dimensions. 

We are interested in potential instabilities when k, k' are near the Fermi surface. 
To be more specific, let k = kp'(/i) be right on the Fermi surface, i.e. ^k = 0, and k' 
be slightly off the surface, with some small ^k', so that A^ = ^k'- See Figure HI L^^^'^ 
is now a function only of the angular directions 6, 9' of k, k', A^, yU and T. As usual, 
the 9 = 7r/4 direction will be referred to as nodal, whereas the 6* = as anti-nodal. 
The chemical potential will be related to the hole doping h by the formula: 

f d'^k 1 

This is clearly an approximation since there are self-energy corrections that modify 
,^k which we ignore; however, since the experimentally measured Fermi surfaces can 
be fit to a ^k of a tight-binding form, we do not expect these corrections to drastically 
affect the main features of our results. Henceforth, we express various properties in 
terms of h defined by the above equation in the limit T — )■ 0. 

Let us first illustrate our findings at the fixed doping h = 0.21. In Figure [5] 
we plot L*^^^) for both k, k' in the anti-nodal direction as a function of A^ at the 
low temperature T = 0.001 (compared to the bandwidth). One sees that for small 
enough A^, i.e. close enough to the Fermi surface, L^'^^'^ becomes positive. As in the 
case of the BCS instability reviewed above, let us identify the gap Agap = A^ with 
the location of the pole in V^g^ in the anti-nodal direction, i.e. the solution to the 
equation: 

— = L('^^)(Agap,/i,T) (18) 
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FIG. 4: Geometry near the Fermi surface in the first quadrant of the Brihouin zone. 

where it is impHcit that 9 = 9' = 0. For Qsc = 1, appropriate to an unscreened (bare) 
g on the order of 10, one sees from Figure [5] that Ag^p ~ 0.003. 
Two other features are apparent from Figure O 

(i) There is a solution to eqn. (TTS!) for arbitrarily large gsc since L^'^^^ passes through 
zero. Furthermore, as gsc — > oo, Agap saturates to where L^'^^'^ crosses the real axis, 
in this case Agap ~ .003. 

(ii) Since L'^'^^^ has a maximum, there are only solutions to eqn. (IT8l) for gsc larger 
than a minimum value, in this case gsc greater than approximately 1/7. Since no 
such threshold was found in 9j, one should conclude that the pairing mechanism 
presented here is essentially different. 

The solutions Agap to eqn. ( fTSll depend on the temperature. However since we 
view T as an infra-red cut-off, one should think of this in terms of the renormalization 
group. Namely, g can be made to depend on T in such a manner as to keep the 
solution Agap fixed. One finds that g increases with decreasing T. For instance, at 
h = 0.15, if gsc = 1 at To = 10^'^, then gsc ~ 2.3 for Tq = 10~^. Henceforth we fix 
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the reference temperature Tq = 0.001, keeping gsc = 1- 
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FIG. 5: The loop integral L'"^^' in the anti-nodal direction as a function of A^ for hole 
doping h = 0.21 and temperature T = 0.001. 

In order to understand how Agap depends on doping, in Figure [6] we plot L^^^'^ 
for h = 0.05,0.15 and 0.24. One sees that as h is decreased, L*^^^'^ crosses zero at 
smaller values of A^, signifying a smaller Agap. Finally, when h is large enough, 
]^iex) jg nowhere positive, signifying no solution to eqn. (ITS]) . In Figure [7] is shown 
the anti-nodal solution Agap as a function of doping h. It is important to note that 
reasonable values for Agap were obtained without introducing any explicit cut-off. 

Numerically, we found no solutions to the eqn. flTSjl in the nodal direction. In fact, 
solutions only exist in a narrow direction around the anti-nodal region. Since the 
anti-nodal direction is along a 1-dimensional chain of lattice sites, this suggests that 
the effect we are describing is closely tied to the properties of the 1-dimensional chains 
studied in the appendix. Although this may help to explain the d-wave nature of the 
gap in the cuprates, it by no means establishes it. In order to study this further, one 
needs a BCS-like gap equation built upon the effective interaction Vgg , or equations 

14 
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h = 0.24 



FIG. 6: The loop integral L^'^^' in the anti-nodal direction as a function of A^ for hole 
doping h = 0.05,0.15 and 0.24. 

analagous to those in 9|, which is beyond the scope of this paper. 

A non-zero t' has a strong influence on the above results. In order to obtain a 
comparable value of the gap, one needs a significantly higher gfsc- For instance, to 
obtain a gap Agap = 0.01 at doping h = 0.15, one needs (^sc ~ 10. 

We wish to point out that the function L^^^^ off-shell, i.e. with A^ equal to an 
arbitrary frequency u unrelated to q appears in the RPA expression for the dielectric 
response function £rpa(w, q) = 1 — gL^^^\A^ = w, q) 13|. Here, u is the frequency of 
an external probe, such as an electric field. Plasmons, i.e. quantized electric charge 
fluctuations, are manifested as delta-function peaks in —lm.{l/efiPA) as a function of 
u. The pairing mechanism studied here thus appears closest to the idea of plasmon 
mediated superconductivity p, [l6|. However our analysis differs in important ways: 
no low energy plasmons were postulated, and the properties of the on-shell L^'^^' are 
very different from those of £:rpa- 
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FIG. 7: Agap in the anti-nodal direction as a function of hole doping h. 
IV. CONCLUSIONS 

In summary, we identified the possibility of Cooper pairing instabilities which 
arise as poles in certain Green functions in the two dimensional Hubbard model, 
in a manner analagous to the BCS theory. Reasonable magnitudes for the gap in 
the anti-nodal direction were found, Agap/t ~ 0.01, without introducing an exphcit 
cut-off. 

In the BCS theory, the analogous poles are 'resolved' by the BCS gap equation. 
The latter has not been developed in the present case, thus this should be investigated 
as a next step in exploring the consequences of the pairing mechanism identified in 
this work. Since it is not clear that the usual BCS gap equation with the effective 
pair potential studied in this paper is valid, this has not been pursued here. Certainly 
the very existence of these poles imply that the effective pairing potential can change 
sign thereby becoming attractive, thus non-zero solutions to the proper gap equation 
should exist. 
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VI. APPENDIX: CHAINS: THE LOOP INTEGRALS IN ONE DIMEN- 
SION. 

In this appendix we study the one-dimensional version of the integrals for L^^^\ 
i.e. eq. f lT^ with J(i^p/(27r)^ — )■ J dp/2TT where p is now a one- dimensional vector, 
and CJk = — 2cosk. In the limit of zero temperature, the occupation number / 
becomes a step function, and the integral can be performed analytically. Namely, 

L(^^) =L+-L_ (19) 

where 



(20) 
where k^(/i) = arccos(— /i/2) is the Fermi momentum. In the region of small q that 
we are interested in, the appropriate branch of the above integrals are given in terms 
of the functions 

,, , 1 / (A£ + 2-2cosq)tan(p/2) -2sinq 
/(p) = arctan — ^ .^ 

7rJA|-8(l-cosq) y JA^ - 8(1 - cosq) 

' (21) 
as follows: 

L+ = /(kp(^)) - /(-k^(/i)), L„ = J(kp(^) - q) - /(-k^(^) - q) (22) 
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In the limit of zero temperature, L^"^^^ is then only a function of q, A^ and the 
chemical potential //, by using the identities 



, /I / WON /2 + /^ , ^n ^ \ ^/o^ V2 + ^ - ^2 - /itan(q/2) 
tan(kH/.)/2) = ^^^, tan((k.(/.)-q)/2) = ^, _ ^ ^ ^, ^ ^tan(q/2) ^^'^ 

Using the above expressions, one can show that L^*^^^ can indeed be positive, which 
is required for pole singularities in V^^ . For instance, for small q'/A^, L^'^^^ has the 
following asymptotic form: 



L(-)(g, A„ /.) ^ ^^^g^ ( 1 + ^ (4 - /^^ - ^1/12) ) (24) 



Whereas q and A^ were treated as independent in the above formulas, physically 
we are interested in the situation where q = k' — k and A^ = ,^k' — ■Ck? referred to 
as 'on-sheir above. In this case there is a divergence in L^*^^^ which arises from the 
pole in the integrand when p = k. Since we are interested in k, k' near the Fermi 
surface, let us be more specific and let k = iipif^) be right on the Fermi surface, 
and k' be slightly off of it, for a small ^k'- Here A^ = ^k' and the singularity occurs 
at qp = arccos(— (^k' + Ai)/2) — ki?(/i). For ^k' = 0.01 and /i = —0.618, which 
corresponds to hole doping h = 0.2, the singularity is at q^? = ±0.0053, which is 
shown in Figure [HI Analytically, this divergence arises as arctan(2) in the above 
formulas. 

The p-integral for L^^^'^ should thus be understood as the Cauchy principal value 
(PV), i.e. /dp — )■ hm^^o (/_^ dp + J^ ^^dp). For finite temperature T, one 
can check numerically that this PV integral is well-defined and finite. However it is 
important to note that L^'^^^ continues to be divergent as T — > 0, which the above 
analytic formulas demonstrate. The reason is that in the limit of zero T, / is a step 
function and the J^ ^ dp piece of the PV is absent, leading it to be ill-defined. This 
is interpreted as an infra-red divergence that needs to be regulated by a finite T, 
however small. 
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FIG. 8: The loop integral L^^^> as a function of q for for A^ = 0.01 and hole doping 
h = 0.20. 

Let now turn to the 1-d analogs of the gap Agap defined in section III. In this Id 
case, because of the divergence as T — )■ 0, it makes sense to define Agap as solutions 
to eqn. (fTSl) with T = Agap, i.e. to determine Agap at a temperature comparable to 
it. The result is shown in Figure [9|, where Qsc was taken to be the screened coupling 
for a bare g = 1. In Figure [TD] we plot Agap as a function of the bare coupling 
g; here, contrary to the 2d case, there is no minimal value of g required for the 
existence of solutions. On the other hand the Agap saturates with increasing g as in 
the two-dimensional case. 



[1] V. J. Emery, Theory of high Tc superconductivity in oxides, Phys. Rev. Lett. 59 (1987) 

2794. 
[2] D. J. Scalapino, E. Loh, and J.E. Hirsch, d-wave pairing near a spin-density-wave 

instability, Phys. Rev. B34 (1986) 8190. 



19 



A 



gap 



0.05 



0.04- 



0.03 - 



0.02 



0.01 




(ex) 

FIG. 9: The Agap pole in V^g as a function of hole doping h for bare coupling g = 1. 

Ap 



0.08 h 



0.06 - 



0.04- 



0.02 - 



igap 




FIG. 10: The Agap pole in V^^ as a function of coupling g for hole doping h = 0.15. 

[3] P. Monthoux, A. V. Balatsky and D. Pines, Toward a theory of high-temperature 
superconductivity in the anti-fermomagnetically correlated cuprate oxides, Phys. Rev. 
Lett. 67 (1991) 3448. 



20 



[4] A. V. Chubukov, D. Pines, and J. Schmalian, A Spin Fluctuation Model for d-wave 



Superconductivity, arXivxond- mat/0201 140 



[5] J. Ruvalds, Plasmons and high-temperature superconductivity in alloys of copper ox- 
ides, Phys. Rev. B35 (1987) 8869. 

[6] C. M. Varma, S. Schmitt-Rink, and E. Abrahams, Charge Transfer Excitations and 
Superconductivity in "Ionic" Metals, Solid State Commun. 63 (1987) 681. 

[7] J. Hirsch, S. Tang, E. Loh, and D. Scalapino, Pairing Interaction in Two-Diniensional 
Cu02, Phys. Rev. Lett. 60 (1988) 1668. 

[8] P. W. Anderson, The resonating valence bond state in La2Cu04 and superconduction. 
Science 235 (1987) 1196. 

[9] S. Raghu, S. A. Kivelson and D. J. Scalapino, Superconductivity in the repulsive Hub- 
bard model: An asymptotically exact weak-coupling solution, Phys. Rev. B81 (2010) 
224505. 
[10] W. Kohn and J. M. Luttinger, New Mechanism for Superconductivity, Phys. Rev. 

Lett. 15 (1965) 524. 
[11] P. A. Lee, N. Nagaosa and X.-G. Wen, Doping a Mott Insulator: Physics of high 



temperature superconductivity, Rev. Mod. Phys. 78 (2006) 17 [cond-mat/0410455| . 
[12] J. R. Schrieffer, Theory of Superconductivity, Addison- Wesley, 1964. 
[13] G. D. Mahan, Many-Particle Physics, Plenum Press 1990. 
[14] A. LeClair, Superconductivity in the two-dimensional Hubbard model based on the exact 

pair potential, larXiv: 1008. 51161 
[15] F. H. L. Essler and V. E. Korepin, Complete Solution of the one- dimensional Hubbard 

Model, Phys. Rev. Lett. 67 (1991) 3848. 
[16] G. D. Mahan and Ji-Wei Wu, Plasmons and high-temperature superconductivity, Phys. 

Rev. D 39 (1989) 265. 



21 



